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Abstract 

The production rates of heavy quarkonia in ion-ion collisions provide sensitive probes in the studies of the hot and 
dense matter formed in these collisions at high energies. However, a reference for understanding the behavior in the 
hot medium is necessary; p-A collisions open the possibility to study heavy quarkonia states in a smaller system of 
much lower average temperature. This is an important step in forming a baseline for understanding A-A collisions, as 
well as an investigation into the nature of modifications of the parton distributions in the nucleus. Using data collected 
at the LHC in 2013, we show results on the prompt J/T^ and ^^(28) nuclear modification factors and the double ratio, 
^^(28) divided by J/'T in p-Pb divided by the same in p-p, in p-Pb collisions at 5.02 TeV. The charmonia states were 
reconstructed via the dimuon decay channel and the yield is analyzed differentially in bins of transverse momentum, 
rapidity, and event activity. 
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1. Introduction 

The suppression of quarkonia in heavy nuclear sys¬ 
tems has been a topic of study for four decades. The 
production rates of J/T^ and ^^(28) mesons from heavy 
nuclei have been measured and compared in the high¬ 
est energy fixed-target and collider experiments since 
1990 Q. 
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Figure 1: The self-normalized ratios for the five measurements from 
ATLAS and CMS listed in the figure key vs. the self-normalized ratios 
of forward calorimeters in each respective experiment. The data have 
not been corrected for any centrality bias, see text. 


There are several motivations for study of the pro¬ 
duction of quarkonia in nuclei. While measurements in 
proton-proton collisions have already provided signifi¬ 
cant constraints on the production mechanism for these 
systems, additional information may emerge from their 
interactions in the extended medium available in proton- 
ion collisions due to the much higher effective gluon 
density made available in such systems. Additional 
information on the formation time and dipole-nucleon 
or hadron-nucleon interaction cross sections for these 
mesons can be obtained from understanding the energy 
dependence of the observed absorption. Further, while 
the production rates of heavy quarkonia states in ion- 
ion collisions may provide a sensitive probe of the hot 
and dense system formed in such collisions, a reference 
measurement is needed in the cold medium in order to 
understand the behavior in the hot medium; p-A colli¬ 
sions open the possibility of studying heavy quarkonia 
states in a smaller system of much lower average tem¬ 
perature. This is an important step in forming a baseline 
for understanding ion-ion collisions, as well as an in¬ 
vestigation into the nature of modifications of the parton 
distributions in the nucleus. 

ATLA8 O has recently measured the production of 
the forward-backward ratio of J/T^ production in proton- 
lead collisions and compared J/'T and 'F(28) produc¬ 
tion in these collisions to the rates in p-p collisions m 
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Figure 2: Plots showing the nuclear modification factor RpPb of J/^ meson production in p-Pb collisions compared to those in p-p collisions. The 
upper two plots are for prompt J/T^ mesons while the lower two plots are for non-prompt J/T^ mesons. The left column shows the px dependence 
while the right column shows the y* dependence. 


The latter comparisons are the focus of the present dis¬ 
cussion. 

2. Experimental Methods 

The measurement of the cross sections for J/'F and 
^^(28) in ATLAS was carried out using the dimuon 
decay channel. One or more triggered muons in the 
hardware-based first level of the trigger (LI), and two 
muons found in a full scan of the software-based Event 
Filter (EF), were required to trigger the event. In offline 
reconstruction, unlike-sign dimuon pairs were identified 
as quarkonium candidates. A calculation of the invari¬ 
ant mass and the pseudo-proper lifetime was performed 
for each candidate in the pair range 8.5 < px <30 GeV 
and center of mass rapidity |y*| < 1.94. A simultaneous 
weighted fit to the invariant mass and the pseudo-proper 
lifetime was performed using a fit model for these line- 
shapes. The weight of each event consisted of the in¬ 


verse of the product of the efficiency of the LI trigger 
for each muon in the pair, the efficiency of the dimuon 
trigger at the EF, the reconstruction efficiency for the 
pair, and the experimental acceptance in the kinematics 
of each pair. 

The collision vertex is determined using the non- 
muonic charged particle tracks in each event. and 
^^(28) mesons originating from the decay of b quarks 
can have a detached vertex reflected in a larger value 
of the pseudo-proper lifetime, while directly produced 
mesons do not have a detached vertex. This is taken 
into account in the fit using separate ’’non-prompt de¬ 
cay” and ’’prompt decay” fitting terms. The fit model 
consisted of four signal terms and three background 
terms. Each signal term consisted of the weighted sum 
of a crystal ball function and a Gaussian function in 
mass multiplied by a term depending on pseudo-proper 
time. For the latter term, the prompt quarkonium pro¬ 
duction was modeled by a delta function while the non- 
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prompt quarkonium production was modeled by an ex¬ 
ponential decay. This model was convoluted by a dou¬ 
ble Gaussian function in pseudo-proper time to take 
into account the experimental detached vertex resolu¬ 
tion. The prompt background in invariant mass was ad¬ 
equately represented as a constant and as a delta func¬ 
tion in pseudo-proper time, while the non-prompt back¬ 
grounds were modeled by two exponential functions in 
both mass and pseudo-proper time. Variations on this 
model were used to estimate the systematic uncertain¬ 
ties, which constituted the largest systematic uncertain¬ 
ties of the measurement, ranging from 2% to 10% for 
J/'F and from 5% to 80% for 'F(2S) , depending on the 
kinematic bin. 

Two different datasets were analysed. First, p-p col¬ 
lisions at = 2.76 TeV were analyzed. Second, p-Pb 
collisions at ^/s^ = 5.02 TeV were analyzed. In or¬ 
der to compare the yields in the p-Pb collisions to p-p 
collisions, it was necessary to perform an interpolation 
of the cross sections between ^[s = 2.16 TeV and ^[s = 
1 TeV, also including data at = 8 TeV when it was 
available, to have the equivalent cross section for p-p at 
^[s^ = 5.02 TeV. The systematic uncertainties of this 
interpolation were carefully studied and are included in 
the overall uncertainties as appropriate. 

An issue of high interest is whether there is any mea¬ 
surable centrality dependence of any observable. For 
central collisions the path length through the medium 
is longest, potentially accentuating any medium-related 
modifications. The experimental determination of the 
centrality in these collisions is more ambiguous than 
in ion-ion collisions and therefore a number of central¬ 
ity measures have been explored. An assumption un¬ 
derlying the centrality measures discussed here is that 
the centrality is monotonically related to the energy de¬ 
posited in the forward calorimeter (FCal). For more de¬ 
tails, see Ref m. 

3. Results 

An approach to reducing the model dependence of 
centrality studies is to use the ’self-normalized yields’ 
as recently employed by CMS in studies of the T(IS) 
state Fig.[^shows self-normalized ratios for prompt 
J/'F and 'F(2S), for non-prompt J/'F , for Z bosons, and 
the CMS points for T(IS). As can be seen, for the high¬ 
est Ft values presumably corresponding to most-central 
events, all ATLAS measurements fall below the diago¬ 
nal line, and the 'F(2S) is far below the other points. 
Plots of this type, corrected for centrality bias in two 
different scenarios, are given in (H. 
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Figure 3: Ratio of the yield of ^(2S) to the yield of Z bosons as a 
function of transverse energy in the forward calorimeter FCal. This 
method of normalizing the yield is expected to strongly reduce model 
assumptions about the event centrality measure. The decreasing trend 
seen suggests a centrality dependence of the T^(2S) production, indi¬ 
cating a dependence on the nuclear medium thickness. 

The nuclear modification factor RpPb can be defined 
as the ratio of the yield in p-Pb to the yield in p-p multi¬ 
plied by 208. This quantity is shown in Fig.|^for prompt 
and non-prompt J/T^ mesons. Generally the trend is for 
this ratio to be somewhat greater than unity, i.e., the pro¬ 
duction is enhanced, and there is no strong dependence 
on Pt or y*. Since for most points the deviation from 
unity is similar in magnitude to the systematic uncer¬ 
tainty for each individual point, and since the uncertain¬ 
ties are substantially correlated, it is difficult to make 
strong statements about this enhancement, however, it 
is an intriguing trend that should be validated by other 
experiments and future studies. 

Another approach to reducing model uncertainties 
in centrality assignments is to normalize the quarkonia 
yield to the yield of particles that do not strongly interact 
with the medium, such as the vector bosons. In Fig. 
is seen the result of this approach when employing Z 
bosons as a normalization for ^^(28) mesons. The over¬ 
all trend of this ratio is to decrease with increasing FCal 
Ft. Plots of this kind for prompt and non-prompt J/T^ 
mesons can be found in Ref. m, and for those mesons 
the trend in the ratio is consistent with being indepen¬ 
dent of FCal Ft. 

Finally, the double ratio RpPb('F(2S))/RpPb(J/'F) 
takes the strongest advantage of cancellation of system¬ 
atic uncertainties, while the normalization should pro¬ 
duce unity in the limit of no medium effects. As can 
be seen in Fig. in the left panel, the mean value of 
this ratio as a function of y* is greater than unity, and 
in the right panel it can be seen that with increasing 
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Figure 4: The double ratio RpPb(T^(2S))/RpPb(J/'F) as a function of FCal Ex. In the absence of medium effects, this ratio should be equal to unity. 
A downward trend with increasing FCal Ex is seen, beginning with a value greater than unity. 


FCal Ex the trend is a decreasing one, beginning at a 
value above unity and ending at a value equal to or less 
than unity. Thus, these data appear to be consistent with 
an enhancement for peripheral events in 'F(2S) produc¬ 
tion relative to J/T^ production, with a transition toward 
a possible suppression for more-central events. 

4. Conclusions 

Measurements of the characteristics of J/T^ and 'T(2S) 
meson production in p-Pb collisions at V(snn) have 
been performed. The differential cross sections for pro¬ 
duction of these particles in prompt and non-prompt 
modes have been determined and employed to charac¬ 
terize the dynamics of the interaction with the medium, 
including the forward-backward ratio and the non¬ 
prompt fraction O . Further, differential cross sections 
for production of these mesons in p-p collisions at ^/(s) 
= 2.76 TeV have been measured. An interpolation of 
these cross sections, together with measurements from 
ATLAS for p-p collisions at ^fls) = 7 and 8 TeV, has 
been performed in order to compare the p-Pb colli¬ 
sion results for these mesons with p-p collisions at the 
same energy of ^l{s) = 5.02 TeV. The comparison of 
J/'P and ^^(28) meson production in p-Pb with p-p has 
been performed through calculation of RpPb for prompt 
and non-prompt J/'P and for prompt ^^(28) as a func¬ 
tion of pt and of y*. These values are consistent with 
being greater than unity and approximately constant 
within the uncertainties. The dependence of the yields 
on event centrality has been investigated using several 
different centrality measures, including self-normalized 
ratios, normalization to Z bosons, and a double ratio de¬ 
scribed below. Both prompt and non-prompt J/'P are 
consistent with being independent of event centrality. 


while a pattern of moderately decreasing yield of ^^(28) 
with increasing centrality is apparent. The double ra¬ 
tio RpPb('P(28))/RpPb(J/'P) has been measured, compar¬ 
ing ^^(28) with J/T^ , and this ratio is equal to unity in 
the limit of no medium effects. This double ratio sug¬ 
gests enhancement for more-peripheral events, tending 
toward a possible suppression for more-central events, 
and the shape of the trend is very similar to that of the 
Z-normalized centrality dependence for ^^(28) . These 
measurements represent a substantial step forward in 
understanding the dynamics of J/'P and 'P(28) produc¬ 
tion in systems involving the cold nuclear medium in 
QCD, as well as improving our understanding of cen¬ 
trality measures in p-Pb collisions. 
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